Luminescence from highly excited nanorings: Luttinger liquid description 
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We study theoretically the luminescence from quantum 
dots of a ring geometry. For high excitation intensities, pho- 
toexcited electrons and holes form Fermi seas. Close to the 
emission threshold, the single-particle spectral lines aquire 
weak many-body satellites. However, away from the thresh- 
old, the discrete luminescence spectrum is completely dom- 
inated by many-body transitions. We employ the Luttinger 
liquid approach to exactly calculate the intensities of all many- 
body spectral lines. We find that the transition from single- 
particle to many-body structure of the emission spectrum is 
governed by a single parameter and that the distribution of 
peaks away from the threshold is universal. 

PACS numbers: 71.10.Ca, 71.45.-d, 78.20.Bh, 78.47.+p 



I. INTRODUCTION 

Photolumincscence (PL) from zero-dimensional ob- 
jects (quantum dots) is one of the highlights in physics 
of nanostructures which emerged during the last decade. 
Early papers (see, e.g., Refs. and the review arti- 
cle Ref. ||) reported the PL spectra consisting of "zero- 
width" luminescence lines. High surface density of quan- 
tum dots caused an ambiguity in assigning of these lines. 
In the later studies the amission from a single dot was 
resolved. This progressu has permitted the PL spec- 
troscopy of individual dots with controllable exciton pop- 
ulation determined by thfi-,excitation intensity, and also 
at tunable charge statesJa-El 

At low excitation intensity the number of excitons in 
a dot is either one or zero. Then the emission line 
corresponds to the transition between the lowest size- 
quantization levels in the conduction and the valence 
bands. Upon increasing the excitation intensity, the 
number of excitons in a dot, N, can be larger than 
one. This leads to the features in the PL spectrum 
which must be interpreted in temjs of recombination 
within many-exciton complexes. lTEj Different recombi- 
nation processes within a complex result in a multitude 
of the emission lines. One reason for an emergence of 
additional PL lines is that, for N > 2, the carriers consti- 
tuting the complex occupy higher size-quantization lev- 
els. Another reason, is that the the interactions between 
the strongly confined photoexcited carriers lift the de- 



generacies of the final many-body states. The latter 
mechanism of interaction-induced multiplication of the 
emission lines .was addressed _in_Refs. pl] , p^p7| for the 
situations withEJ and withouttlE3 orbital degeneracy of 
single-particle states. The calculations carried out in 
Refs. 0,qJ|,[L7| predicted the splittings of the PL lines, 
originating from different many-body final states, to be 
of the order of the matrix element of the interaction po- 
tential. The actual positions of the lines predicted by 
these calculations reproduce quite accurately the experi- 
mental PL spectra of Refs. p~T] |l4| (for up to N = 16) and 
of Ref. HI (for 1 < N < 6). 

Within the approaches of Refs. |l^JT^,|l^, each many- 
body line corresponds to recombination of an electron 
and a hole having the same size-quantization quan tum 
numbers. We note here that in the experimenttlrEJ ad- 
ditional emission lines have been observed that were 
identified with the transitions between different size- 
quantization levels. These transitions originate from 
a shake-up effect in a confined electron-hole system. 
Namely, the radiative recombination of an electron-hole 
pair is accompanied by an internal excitation within 
the exciton multiplex. These many-body processes 
are missed within the theoretical approaches of Refs. 
[H] , |l2||i7[ Meanwhile, as we argue below, the PL lines 
originating from the shake-up processes multiply very 
rapidly with increasing N, and for high enough N be- 
come the dominant feature of the PL spectrum. Here we 
develop a theory of the many-body luminescence from a 
quantum dot for the limiting case N ^> 1. 

When a number of photoexcited carriers is large, an 
adequate description of PL from a dot must be devel- 
oped in terms of Fermi seas formed by equal numbers 
(determined by excitation) of electrons and holes in con- 
duction and valence bands, respectively (see inset in Fig. 
|l|). As in the case of a small number of carriers, such 
a description is based on the fact that PL is preceeded 
by a fast nonradiative relaxation of electrons and holes 
into the corresponding ground states .Eirt3 Microscopic 
mechanism of this relaxation is addressed, e.g., in Ref. 

For noninteracting system, the emission lines would 
correspond to transitions between size-quantization lev- 
els in conduction and valence bands which obey the se- 
lection rules. Within this description, the single-particle 
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emission spectrum near the Fermi edge is given by the 
Golden rule, 



of peaks within the neighboring generations are almost 
identical. 



J(w) oc ^2 C n5 w + (Ai + A 2 )n 



(1) 



where Ai and A2 are level spacings for electrons and 
holes; C n are the oscillator strength which depend on n 
only weakly (w < is measured from the Fermi edge). 

As discussed above, the many-body transitions, result- 
ing from interactions of carriers in a dot, change qualita- 
tively the form of the spectrum. Analogously to recom- 
bination within an exciton multiplex, here a removal of 
an e-h pair shakes up the respective Fermi seas by caus- 
ing them to emit Fermi sea excitations. Since in a finite 
system, the energies of excitations are quantized, such a 
shake-up would lead to the spectrum of a form 

I(lu) cx C »™5 (w + mAi + nA 2 j , (2) 

ran 

rather than Eq. (|l|). Here Ai and A 2 are the level spac- 
ings renormalized by interactions. All the information 
about many-body correlations in the system is encoded 
in the oscillator strengths C mn . As we will see below, 
Cram being governed by interactions, are strong func- 
tions of to and n. 

The goal of the present paper is to demonstrate that 
the oscillator strengths C mn can be evaluated analytically 
for a dot o£ a ring geometry. Such dots have been recently 
fabricatedQ and their emission spectra (including many- 
body effects) were Sjtudifid for low excitation-intcnsitics 
both experimentallyli3i23 and theoretically.EiTc3 Our ap- 
proach is valid when the number of carriers, N, is large 
enough (the accuracy is 1/N). However, this "asymp- 
totic" consideration allows us to establish the universal 
properties of the many-body spectrum away from the 
threshold. 

For a ring-shaped dot, the electron and hole Fermi seas 
represent one-dimensional (ID) systems. This allows us 
to use the finite-size Luttinger-liquid descriptionEZl for 
calculation of the emission spectrum. Note that the Lut- 
tinger liquid model was employed earlier for calculations 
of the Fermi-edge optical properties of infinite ID sys- 
tems (with and without defects) in Refs. |28|-|32]. 

We show that due to a finite size of the system, the 
structure of the emission spectrum is different in low- 
and high-w domains. Namely, for ^tln|cj/(Ai + Ai) <C 1 
(low frequencies) the spectrum is dominated by single- 
particle peaks; many-body satellites have a relative mag- 
nitude ~ fi, where /1 C 1 is the dimensionless inter- 
action strength (Luttinger liquid parameter). For high 
frequencies (i.e., /Ltln|u;/(Ai + Ai) ^> 1), the many- 
body peaks completely dominate the spectrum; roughly 
speaking, in the high-w domain, the oscillator strengths 
of single-particle peaks are evenly distributed among the 
multitude of many-body peaks. Furthermore, the peaks 
cluster into groups (generations), so that the patterns 



II. EMISSION FROM A LUTTINGER LIQUID 
RING 

We start with— .the two-component Luttinger liquid 
model on a ringpa with Hamiltonian Hi + H2 + Hi nt , 
where Hj describe noninteracting electrons (j = 1) and 
holes (J = 2) with linearized dispersions (the slopes are 
determined by the Fermi velocities Vj); Hi nt describes 
the interactions between carriers via screened potential 
U (x). The e-h recombination rate is given by the Golden 
rule 

W(u J ) = T Y,\(f\m\ 2 5(E l -E f -c) 
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i\T\t)T{Q)\i), 



(3) 



where Ei and Ef are the energies of initial (ground) and 
final (with e-h pair removed) states, and 



T = T + + T_ T ± 



d I dxip2^(x)ipi±(x) (4) 
'0 



is the dipole transition operator. Here ipi± are annihila- 
tion operators for left (— ) and right (+) moving carriers, 
d is the interband dipole matrix element, and L is the 
ring circumference. Note that recombination occurs be- 
tween left (right) electrons and right (left) holes. The 
recombination rate is then expressed via a four-particle 
Green function, 



W(uj) =d 2 dx 



dte 



Jt [D + (x,t) + D-(x,t)] 

= d 2 [D + (uj) + D_(cu)], 
D ± (x,t) = (4 T (x,t)4 ± (x,t)i; 1± (0)^ 2T (0)). 



(5) 



In order to svaktate D±(x,t) for a two-component Lut- 
tinger_liquidJ23i23 we use the bosonization technique on 
a ringE3 (see Appendix A). The final result reads 



D±{x,t) 
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(6) 



where Sp = n(v\ -\-v%)jL is the energy shift (to be ab- 
sorbed into the frequency) due the change in the parity 
of electron and hole numbers, and e is a cutoff. The coor- 
dinate dependence of D a (x, t) is determined by (a = ±) 
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The explicit expressions for renormalized Fermi velocities 
Vj and interaction-induced exponents Hj are given in Ap- 
pendix Correspondingly, the level spacings are now 
Aj = 2irvj/L. The interaction strength is characterized 
by the ratio u/vj, with 



J- [ dxUix) 

hit J 



(8) 



being the Fourier of screened potential; this ratio repre- 
sents the average (screened) interaction in units of the 
(bare) level spacing near the Fermi energy. For weak 
interactions, u/vj <C 1, we have (see Appendix [a]) 



Hi si {u/A Vj f, Aj ~ A,(l + u/2 Vj ). 



(9) 



The correlator D a (x,t) is periodic in variables Zj a . In 
order to carry out the integration in Eq. (^) we first per- 
form the Fourier expansion of [/ Q (£j Q )l as 



/±0j±) 
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B{n + v,l-v)e ±2mz i ±IL , 



(10) 



where B(x,y) is the Beta- function. Substituting this ex- 
pansion into Eq. (^), and then into Eq. (^), we arrive at 
the expected form Eq. (Q) of the emission spectrum with 
the coefficients C mn cast in the following closed form (see 
appendix H) . 



* dfadfodfo e^+^e-sWi+'feHm+n) 
— e^ 1 ] I 1 — e 1 ^ 2 



(27T) 



e 2 



4>3(m-n) 



(11) 



Note, that the sum in Eq. (|2j) is constrained by the se- 
lection rule that m and n are of the same parity, i.e., the 
combinations 



N = (m + n) /2, M = (m - n) /2, 



(12) 



which enter into the rhs of Eq. Qllj), are integers. This is 
the result of the linear dispersion of electrons and holes 
near the Fermi levels. 

Formula ( O ) for the oscillator strengths is the main 
result of this section. It is easy to see that it correctly re- 
produces the non-interacting limit. Indeed, upon setting 
= 0, the integral over (j>3 yields C mn = S mn . An- 
other important limiting case m, n 3> 1 corresponds to 
the transitions well away from the Fermi ed ge. In this 
case, the main contribution to the integral (O) comes 
from the domain <p% + </> 2 ~ (m + n) _1 <C 1. Within this 
domain, one can neglect the difference between <pi and 
—<f>2 in the last two factors in the denominator. Then the 
integrals over (j)i,4>2 factorize, yielding 



C mn =S^ +1 ± MN ll +f t±K(M), 



r(i + / i 1 )r(i + Al2 )[r(A^ + i)] 5 



(13) 



with 



J -TV 27T (l 
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e 2(Mi+W)(_ 1 )M r ( 1 _ fJLl 

r(l-M- Ml )r(l + Af-^ 2 ) ' 



(14) 



where T(x) is the Gamma- function. It can be seen from 
Eq. jl4| ) that, for a given N, the oscillator strengths, 
C mn , fall off as C mn oc |Af | Ml+Al2_1 with increasing 
\M\ — i|m — n\. This slow power-law decay reveals 
strong correlations within electron-hole system on a ring. 
Finally, using the large x asymptotics of r(x), we ob- 
tain the expression for the oscillator strengths valid for 
N, \M\ » 1, 



e 2M r(l — ji) sm-Kp, 



N"\M\"-\ (15) 



r(i + Mi)r(i + M2 ) 7T 

where fi = fi\ + fi 2 , and fi = \[i + - /i 2 )sgnM. 



III. MANY-BODY STRUCTURE OF THE 
EMISSION SPECTRUM 



The general expression (|llj) determines the heights of 
the emission peaks, while the order of the peaks with 
different {m, n} is governed by the 5-functions in Eq. 
(^|), which ensure the energy conservation. Therefore, 
this order depends crucially on the relation between Ai 
and A 2 . Moreover, a commensurability between Ai and 
A 2 leads to accidental degeneracies in the positions of the 
emission lines. However, in order to establish the general 
properties of the spectrum, it is instructive to consider 
first two particular cases of commensurate Ai and A 2 . 

We start with the symmetric case A, = A/2 (and, 
hence, fa = ;u/2). The peak positions, as determined by 
Eq. (||) , coincide with those for single-particle transitions, 
| a; | = NA. The corresponding oscillator strengths can be 
straightforwardly evaluated from Eq. (O) as 



cn — ^ Cn+m,n-m — 

M 



e fJ. e -iN4> 



-i 2 



2tt 



1 



1+2/i 



(16) 



For N 3> 1, the denominator of the integrand can be 
expanded, yielding 



CN ^ (eN) 



2a _ 



1 



2/i 



(17) 



Note that single-particle oscillator strengths correspond 
to cat = 1. We thus conclude that interactions af- 
fect strongly the peak heights for |( J j/A| 2 ' 1 3> 1, i.e., 
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in the high frequency domain. In fact, even for an ar- 
bitrary relation between Ai and A2, the crossover be- 
tween "single-particle" and "many-body" domains of the 
spectrum is governed by the dimcnsionlcss parameter 
/xln|w/(Ai + A 2 )|. 

Indeed, consider now the case Ai = 3A2 (and thus 
iii ~ 9/ii) which renders a spectrum richer than ([!]) and 
(|l7|). As follows from Eq. (||), the spectral positions of 
the peaks are given by |w|/(Ai + A 2 ) = 1/2, where I is 
an integer. The corresponding oscillator strengths can 
be evaluated explicitly from Eq. ([ll]) (see appendix |^) . 
The final result reads 



D±(w) = — 



A 



2/' 



I . 



1 + \~ 



2u I 
A I 



■ 8{uj + AO 



1 — 1 — 1 



cu + A(l + 1/2) 



(18) 



with fi = fii + [j,2 — 10/Ki and A = Ai + A2 = 4Ai. 

The above result illustrates how the structure of the 
spectrum evolves as the frequency departs from the Fermi 



For yt/. In I — ■ I <C 1, each single-particle peak, \u\ 



edge. 

I A acquires a weak many-body satellite at |w| = (l+-k)A. 
In the opposite limit, /iln| # | ^> 1, the oscillator strength 
of an "integer" peak gets equally redistributed between 
the components of the doublet. The crossover frequency, 
ft, separating the "single-particle" and the developed 
many-body domains of the spectrum is determined by the 
condition ln| -2. | ~ fjT 1 . The spectrum ( |l8| ) is schemati- 
cally depicted in Fig. |l|. 

Let us turn to the structure of the spectrum in the 
general case of incommensurate Ai and A2. We start 
from the observation that the peak positions can be clas- 
sified by "generations" . Namely, once a peak {to, 0} 
(or {0, n}) emerges at u> = to m — — mA\ (or uj = 
uj n = — nA 2 ), it is followed by next generations of peaks 
win' = w m - fc(Ai + A 2 ) or w^' — uj n ~ k(Ai + A 2 ) 
repeating with a period A = Ai + A 2 . Thus, for a crude 
description of the spectrum away from the Fermi edge it 
is convenient to divide the frequency region to < into 
the intervals of width A. 

The number of peaks within the spectral interval 
{— \uj\, —\w\ — A} is the number of integers satisfying the 
conditions \w\ < mAi + ?iA 2 < |cj| + A. This number is 
equal to 



AC = 



2AxA 2 ' 



(19) 



where we assumed \w\ 3> A and took into account the 
parity restriction. From Eq. ( |l9| ) we find the peak density 
9u = AC/A = |w|/2A 1 A 2 . It also follows from @ that 
A/" w _^ — AC = A 2 /2AiA 2 generations start within each 
interval. Since the heights of consecutive peaks within 



the interval A vary non-monotonically, it is natural to 
characterize these heights by the distribution function 

1 f°° - 
F(C) = - — / dmdnS(uj + toAi + nA 2 ) S(C mn - C), 

(20) 



where C mn is given by Eq. (15). Here we made use 
of the fact that AC ^ 1 by treating to and n as con- 
tinuous variables. The prefactor in Eq. (|20| ) ensures 
the normalization (J Q dCF(C) = l). It is easy to see 
that F (C) is nonzero in the interval between C m in = 

min{2^i j2 |% 1 | 1 2m } and C max = 2| ^ | M max{^i :2 } (we 
omit the overall factor e 2AI ). Within this wide interval, 



F(C) falls off as (C /C) 2+M , where 



C = i-i 



2/^-1 



(21) 



is the typical value of the oscillator strength. On the other 
hand, the average oscillator strength, which can be easily 
calculated from Eq. (p0|), is equal to C = ^~ 1 Cq S> Co- 
The distribution function F(C) is schematically depicted 
in Fig. S. The fact that C decreases with \uj\ can be 
understood in the following way. As it is seen from Eq. 
(|i"7|), in the symmetric case, with only a single peak per 
interval A, the peak heights increase with \oj\ as |^| 2m . In 
the general case, this spectral intensity gets redistributed 
between AC different peaks. Thus, 



C~AC 



A 



2/i 



OC \UJ 



2/i-l 



(22) 



IV. CONCLUSIONS 

In the present paper we derived the emission spectrum 
from a highly excited ring-shaped quantum dot. In this 
system electron-electron, hole-hole and electron-hole in- 
teractions relax the momentum conservation leading to 
a multitude of discrete emission lines. Luttinger liquid 
model employed in our calculation allows to evaluate the 
overlap integrals between the correlated initial and final 
many-body states. These overlap integrals determine the 
intensity of the corresponding spectral lines. 

The theoretical value of the dimensionless interaction 
parameter /1 is determined by the ratio of screened in- 
teraction U to the level spacings Ai and A 2 at the cor- 
responding Fermi levels. Both quantities depend on the 
number of excited carriers, N , which in turn is deter- 
mined by the excitation intensity. This, and the sensitiv- 
ity of the screening to the details of experimental setup, 
lead to a common ambiguity in the theoretical determi- 
nation of /i. For example, in quantum wires, th e Vfil ue of 
jj, measured in resonant tunneling experiments,! 3 ^r 8 ! was 
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significantly larger then theoretical estimates. Concern- 
ing the estimates for Ai and A2, in the experimental 
paper Ref. || on luminescence from ring-shape dots, the 
total energy separation A between the lowest level was 
approximately 5 meV. This value comes almost exclu- 
sively from the conduction band, due to the large ratio 
of the electron and hole effective masses. Both Ai and 
A2 increase linearly with increasing N. This implies that 
the shake-up processes within the hole system are exper- 
imentally much more relevant than those for electrons. 

Note finally, that for emission from a finite electron- 
hole ID system considered here, the physics underlying 
the interaction-induced multiplication of the number of 
lines with departure from the Fermi level is analogous to 
that for tunneling into a disordered quantum dotrJ 
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APPENDIX A: 

Here we outline the calculation of the Green function 
(^J) using a bosonisation scheme for the multicomponent 
Luttinger liquid on a ring.LJ The right/left fermion fields 
are presented as 



ip ja (x) = (2 7r e)- 1 /2 e ^ Q (^)+^/i ; 



(AI) 



where right/left (a — ±) bosonic fields (pj a (x) are related 
to the corresponding densities as Pj a {x) — ^ d{ Po»( x ) 
(here e is a cutoff). The bosonic field has a decompo- 
sition 



tpj a (x) = (p° ja + aN ja 2irx/L + ipj a (x), 



(A2) 



where the number operator Nj a and its conjugate (p° Q 
satisfy the commutation relations 



(A3) 



and the periodic fields ipj a (x) — <fij a (% + L) have the 
usual form. 



<p ja (x) = E^ a )\/i^j e " l9|e/2 ( ei9Xa « + 



(A4) 



with a q j and a/ q . satisfying standard boson commutation 
relations \9(x) is the step function]. The boundary con- 
dition for the fermion fields, i/jj a (x+L) = (—l) Nj ^pj a (x), 
depends on the parity of the number of particles, Nj = 
2Nj a . The Hamiltonian H — Hq + Hi n t is quadratic in 
boson fields: 



/ dx 


dipj a (x) 


la 


dx 



and 



Hint= -2^l dx L 



dy 



E 



a d(pj a (x) 



E 

L 



2tt dx 

d<pip(y) 



(A5) 



Uji(x-y) 



27r dy 



(A6) 



wher e Uj i(x) is the screened potential. Using Eqs. ( [A2| ) 
and (A4), and after separating out the zero-mode part 
of the Hamiltonian, H°, from the bosonic part, H, the 
total Hamiltonian H = H° + H can be written as 



H = 7 E NjaivjSjlS, 



jla/3 



Oa/3 ■ 



Ujl_ 

2 



N, 



iff 



E 



e~^\q\ 



VjSjiaLaqj 



+ M(a+.+a_ gi )(a t _ 



ql 



I'll, 



(A7) 



where Uji = ir 1 J dxUji(x). 

In order to calculate correlation functions, the Hamil- 



tonian H, corresponding to the second term of (A7), has 
to be brought to the canonical form. This is done in two 
steps. First, we perform a two-component Bogolubov's 
transformation in order to eliminate the cross-terms with 
opposite momenta, 



a u = ^2( X M b ii + Y ji b - q i)> 
1 

We then obtain 
H = £ e-^|< ( X t _ y t).^ (x „ Y)inhqr 

where the matrices X and Y must satisfy 
J2(xi + Yi). l (u lm +v l 6 lm )(X + Y) rnn 

lm 

= U x '- Y ') 3 M x - Y ) ln - 



(A8) 



(A9) 



(AlO) 



Second, we diagonalize the Hamiltonian ( A9) by first pre- 
senting the matrices X and Y as 
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A = coshAO, y = sinhAO, 



(All) 



where Xji = Xjdji is diagonal matrix of Bogolubov's an- 
gles Xj and O is an orthogonal matrix, and then by in- 
troducing new boson operators c q j — J2i^ji^qi- The 
Hamiltonian H then takes the form 



H = Y,e-^\q\v jC \ jCqh 
qj 



(A12) 



-2A, 

old and new boson operators are related as 



a q j = cosh Xj c q j + sinh Xj c_ 



'i j ' 



(A13) 



Using the decomposition (All), Eq. ( A10 ) takes the form 
OAO = 0, where O is the transposed matrix, and the 
matrix A is given by 

An = Uj ie^ +Xt + Sjivj (e 2 ^ - e" 2 ^) . (A14) 

The Bogolubov's angles Xj are found from the condi- 
tion that all the eigenvalues of Aji vanish. In the two- 
component case, this yields 



-2Ai _ 



Q 




= Q/e- 2X \ 



V\V 2 



(A15) 



The Luttinger liquid interaction parameter is given by 
fij = sinh 2 Xj. In the case of weak interactions, Uji/vj <C 
1, we have Xj ~ —Ujj/Avj so that fij ~ A 2 ~ (ujj/Avj) 2 

and Aj ~ Aj(l + Ujj/2vj) 



With the Hamiltonian (A12), the time-dependence of 
new operato rs is standard, c q j(t) — e_^ j ' ? '*c gj -. Using 



the relation (A13), the periodic fields (A4) take the form 



<p ja {x,t) = 




q\e/2 



cosh A, 



</./ 



'(— qa) sinh Xj 

a -iqx+iVj\q\t J 

qj 



(A16) 



The time-dependence of zero-modes is governed by the 
zero- mode part (first term) of the Hamiltonian ( |A7| ) . The 
time-dependent bosonic field is finally obtained as 

cpj a (x,t) — ip° a + aNj a 2ir(x - aVjt)/L 

-^2ujiNifiirt/L + (p ja (x,t). (A17) 
10 

We are now in position to calculate the Green func- 
tions. For this, we separate out annihilation and creation 
parts of the periodic field ( Al6| ), (fj a (x,t) = (pj a (x,t) + 



(p~j a (x,i), which satisfy the following commutation rela- 
tions 



with 



In f a (z ja - z! ja )f- a {zj- a - Zj_ a ) , (A18) 

x — aiijt. Then we present the fermion oper- 
ator (Al) in the normal-ordered form, 

i)ja(x,t) = 1p° a (x,t)lpj a (x,t), 

^^iaN ja 2nz ja /L-i^2 ifj ujiNip nt/L 

$ ja (x,t) = L-^hne/Ly^^U^e^I*^, (A19) 

where we again s epar ated out zero-mode and peri odic 
parts. Usi ng Eq. ( Al9| ) together with commutators (A3) 
and ( [A18D , the Green function @ can be straightfor- 
wardly calculated as 



D a (x,t) 



2 ?re \ 2(M2+/J2) e -it8 P -itS u 



L 



fa(Zla) 



1+Mi 



L 2 

f-a(zi-a) 



f-a(z2,-a) 



1 + M2 



fa(z 2a ) 



(A20) 



where 5p = ir(vi + V2)/L and S u — n(uii + U22 + 2ui2)/2 
are the energy shifts due the changes in the parity of 
electron and hole numbers and in the Coulomb energy, 
caused by a removal of an e-h pair. We assume that the 
screened interaction is the same for electrons and holes, 
Mil = U22 = — u\2 = u, so that S u = 0. Then, after- 
absorbing the factor 2n/L into e, we arrive at Eq. (^J). 

Note finally that the above calculation is easily general- 
ized if the ring is penetrated by a magnetic flux <f>. In this 
case, the electron and hole number operators should be 
shifted by flux-dependent constants, N\ a — > N la +ot(j)/ fa 
and N201 ~ ot<f>/ '4>o, where </>o is the flux quantum. 

This results in a replacement Sp — > Sp(l — 2a4>/4>o) in 
Eq. (§). 



APPENDIX B: 



Substituting the Fourier expansion 



bu{n) = B(n + v, 1 - v), 



into Eq. (p), D a (uj) takes the form 



(Bl) 



M 

xb^(n' 2 )A ± (Lu,{n}), (B2) 



G 



with 

A±(uj,{n})-- 



1 

1? 



dt I dxexp 
o 



-iojt 



2?r 
2tt 



' — ("i ^1 : - n[z lT - n 2 z 2T + n' 2 z 2 ±^j 



— 7l\ ,71-2— TlU ^ 



(n 2 + n' 2 ) , (B3) 



where we absorbed the parity shift Sp into u. The Kro- 
niker delta and the delta-function reflect the conservation 
of momentum and energy, respectively. Thus, we obtain 



D a (u) = ^E + ^ lTn + ^ n )- 



(B4) 



with 



C mn = e 2 ^+^> J2 [("* + n)/2- l}b 1+ ^ 2 (n - I) 
i 

x6 w [(m-n)/2 + l]6 Pa (l). (B5) 

Finally, using t he in tegral representation for the Beta- 
function in Eq. (Bl) we arrive at Eq. (|ll|). The sum in 
Eq. (B4) is constrained by the selection rule that m and 
n are of the same parity, as can be seen from Eq. (B2). 
From Eq. (B4), the emission spectrum (0) follows. 



APPENDIX C: 



Here we consider the case when the level spacings in 
the conduction and valence bands are commensurate: 
A1/ A2 = p/q, where p and q are integers. Then Eq. 
(B4) takes the form 



D a (w)= — 2J C mn S (uj + A 



~ mp + nq 



p + q 



2tt 



Y,C k 5(u + Ak/Q), 



(CI) 



where Q = p + q, P = p — q, A = Ai+A2, and 



Ck= E h 



Uk,MP+NQCN+M,N-M 

mn MN 



E^ 



AfiV 

Using the relation 



(C2) 



1 

E'w = q E e 



the oscillator strengths can be presented as 



(C3) 



Q-i 



^ = AE e ^ fc/Q /'( fc )' 



(C4) 



with 



/i (fc) = y e ^MP/Q C , (C5) 



Using integral representation (11), the sum over M can 
be explicitly performed. For k/Q — \u>\/A 3> 1, the 
resulting expression for coefficients fi takes the form 

F 2(^i+/i 2 ) p -i(<t>i+4> 2 )k/Q 

fl(k)= 



, (2tt) 2 (— i^i)^^ 1 (-i0 2 ) 1+Ma 
1 

(1 - s; - i Sl [<j) 2 + (4>i + 4> 2 )P/Q}Y 
1 

(l - af - iaftfa - {fa + (MP/Q]) 1 



-, (C6) 



where si = e l27Tlp /Q . The ^-dependence of /z(fc) is deter- 
mined by the relative magnitude of Q/k and 1 — si\: 

2( J U1+ J U2) 

(C7) 



/«(*) 

for fc/Q < |1 - s;|~\ and 



cfc 



e 2 fc/Q 

l-Si 



e 2 k/Q 
1-st 



(C8) 



for fe/Q 3> |1 — s/| 1 , with the two estimates matching 
at fc/Q~ ll-Sil" 1 ^. 

In the case A1/A2 = 3, corresponding to P = 2 and 
Q = 4 so that s; = (—1)', the coefficients fi take two 
different values depending on the parity of I, 



feven (&) 

fodd(k) : 
yielding 

Cfc = 



2(^1 +/i 2 ) 



A»l+A»2 



cu; 



2(fil+/i 2 ) 



2A 



M1+M2 



2p 



l + (-l) fc l + e l7rfe/2 |^| " 



(C9) 
(CIO) 

(en) 



with n = /ii + /i2- Obviously, Cfc = for A: odd. For A; 
even, we have 



Gy 



a" 



2(ii 1 I 2u> I ~M 
I A I 



G 



4i+2 



2(11 _ I 2a; I -A" 
I A I 



(C12) 
(C13) 



leading to Eq. ( |18[ ) . 

*On leave from Vanderbilt University 
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FIG. 1. Emission spectrum for Ai = 3A2. In the 
low-frequency domain, the single-particle peaks acquire weak 
many-body satellites; in the high-frequency domain, the 
heights of the single-particle (|w|/A = V) and many-body 
(|w|/A = I + 1/2) peaks are close to each other. Inset: Sin- 
gle-particle spectra of electrons and hole in conduction and 
valence bands, respectively; u t h is the energy distance be- 
tween the corresponding Fermi levels. 



FIG. 2. The distribution function Eq. ( pfj| ) of the peak 
heights within the interval A is plotted schematically ver- 
sus x — C/Cq. The minimal value of x is x m i n ~ 1, while 
x m ax ~ l^/Aj M ^> 1. The point x = corresponds to 
the average oscillator strength. 
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Fig. 2 
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